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INTRODUCTION 


Radar  detection  models  vary  from  simple  models  (page  20,  reference  1)  with  minimal 
refinements  to  more  sophisticated  models  (reference  2)  with  many  effects  consldezed, 
e.g.,  atmospheric,  earth's  surface,  target,  receiver,  operator,  clutter.  It  Is  reasonable 
to  e^ct  that  with  more  refinements,  a model  can  more  closely  predict  the  real  slgnal-to- 
noise  (S/N)  ratio  under  the  assumed  conditions.  It  Is  felt  that  the  primary  difference 
between  model  predictions  and  fleet  results  lies  In  the  modelli^  accomplished  between 
the  determination  of  the  S/N  ratio  and  the  resultant  probability  of  detection. 

This  paper  does  not  delve  Into  the  probablUty  of  detection.  Instead,  the  paper  pre- 
sents a radar  detection  model  that  calculates  the  S/N  ratio  as  a Amctlon  of  target  range. 
The  user  must  decide  what  techniques  to  apply  when  converting  the  S/N  ratio  to  the  pro- 
bability of  detection. 

The  model  Is  a revision  of  the  radar  detection  model  from  the  "Countering  the 
Anti -Ship  Missile  (CAM)  Study, " reference  3.  There  are  three  additional  capabilities: 

e Self -screening  Jammers 
e Ascending/descending  target  flight  profiles 
e Simultaneous  main-lobe  and  side-lobe  Jamming. 

The  CAM  model  was  further  refined  by  improving  the  methods  used  to  calculate  the 
multipath  effect  and  atmospheric  attenuation. 
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RADAR  DETECTION  MODEL 


GENERAL 

A macro  flow  diagram  of  the  revised  radar  detection  model  is  presented  in  figure  1. 
The  required  inputs  to  the  model  are  briefly  discussed  in  table  1.  The  detection  model, 
in  FORTRAN  IV  language,  is  located  in  appendix  D with  user  instnictions  in  appendix  A. 
A standoff  Jamming  example  (appendix  B)  and  a self-screenlqg  jamming  example  (appen- 
dix C)  are  provided  to  assist  users  in  running  the  program. 

MODEL  STRENGTHS 

The  model  determines  S/N  ratio  as  a function  of  target  range  In  a clear  or  noise 
jamming  environment.  The  noise  jamming  may  be  created  by  a standoff  or  self-screening 
jammer.  The  threats  may  have  an  ascending,  level,  or  descending  flight  profile.  The 
effects  of  atmospheric  attenuation,  weather  clutter,  surface  clutter,  sea  state,  and 
multipath  may  be  considered. 

In  the  multipath  calculations,  the  model  has  a significantly  improved  method  of  cal- 
culating the  reflected  ray  path  length. 

The  a and  /3  ranges  are  commonly  used  measures  of  radar  performance: 

Q range  — the  detection  range  against  a one  square  meter  target  with 
a one  watt  per  megahertz  self -screening  jammer. 

^ range  — the  detection  range  against  a one  square  meter  target  in  a 
clear  ECM  environment. 

The  model  can  calculate  the  a and  0 ranges  to  provide  the  analyst  with  a test  that 
the  radar  being  examined  is  being  modeled  correctly  by  comparing  the  a and  ranges 
to  values  obtained  from  other  reliable  sources. 

MODEL  WEAKNESSES 

The  model  calculates  the  S/N  ratio  as  a function  of  target  range.  Then  the  user 
must  convert  the  S/N  ratio  to  probability  of  detection  using  appropriate  techniques. 

The  targets  must  be  opening  or  closing  the  radar's  location  on  a radial  bearing. 

The  model  does  not  allow  for  an  offset  radar  site  that  produces  crossing  targets. 

Refractivity,  the  bending  of  the  radar  waves  by  the  atmosphere,  is  accounted  for  by 
modifying  the  radius  of  the  earth.  The  amount  of  refraction  may  be  varied,  but  special 
effects  - e.g. , surface  ducting,  layering  - may  not  be  considered.  The  amount  of  re- 
fraction is  programmed  using  input  X(15):  the  coefficient  of  refractivity. 
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FIG.  1:  RADAR  DETECTION  FLOW  DIAGRAM 
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X(7)  and  X(8) 
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RADAR  DETECTION  MODEL  INPUTS 
DEFINITION 

multipath  effect  Included  (Isyes;  0=no) 

surface  clutter  Included  (l=yes;  O=no) 

surface  clutter  coefficient  In  decibels.  A ratio  of  1 
square  meter  of  clutter  to  1,000  square  meters  illu- 
minated would  be  entered  as  -30  (see  reference  1, 
pages  527-534) 

weather  clutter  Included  (l=yes;  0=no) 

rainfall  rate  in  millimeters  |)er  hour 

attenuation  included  (l=yes;  0=no) 

X(7>=L^  and  X(8)«L2  of  natural  logarithmic  regression 

L=Lj  + L2  LnR,  where  L is  the  attenuation  loss  in 

decibels  and  R is  target  range  in  nautical  miles  (see 
section  on  attenuation) 

jamming  including  (l=yes;  0=no) 

type  of  jammer  (0=standoff  jammer;  l=self-screening 
jammer) 


X(ll) 

X(12) 

X(13) 

X(14) 


I 


ii  . 


I ! 
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location  of  jamming  (0=side-lobe;  l=main-lobe:  2=both 
main-  and  side-lobes) 

range  from  radar  to  jammers  in  nautical  miles.  Note: 
all  jammers  must  be  at  the  same  range 

noise  jammer  bandwidth  in  megahertz.  Note:  all  jammers 
must  have  identical  bandwidths 

effective  radiated  power  of  the  jammer(s)  in  watts. 

Packed  word  "SSSSMMMM"  where  "SSSS"  refers  to 
side-lobe  jamming  and  "MMMM"  refers  to  main-lobe 
jamming.  Note:  all  jamming  power  in  the  appropriate 
lobe  must  be  totaled  and  entered  as  one  value.  For 
example,  three  100  watt  jammers  in  the  side -lobe  and 
two  100  watt  jammers  in  the  main-lobe  would  be  entered 
as  03000200 


-4 


X(33) 

X(34) 


TABLE  1 (Cont’d) 

DEFINITION 

coefficient  of  refractivity  (1.333  for  the  4/3  earth 
approximation) 

sea  state  (Beaufort  scale) 

target  altitude  at  the  beginning  of  the  run  (in  feet) 

target  altitude  at  the  end  of  the  run  (in  feet) 

target  range  at  the  beginning  of  the  run  (in  nautical  miles) 

target  range  at  the  end  of  the  run  (in  nautical  miles) 

average  target  cross  section  (in  square  meters) 

peak  radar  transmitted  power  (in  kilowatts) 

main-lobe  antenna  gain  (in  decibels) 

side-lobe  antenna  gain  (in  decibels)  (e.g. , side-lobe 
down  25  db  from  the  main-lobe  gain  would  be  entered 
as  25) 

radar  frequency  (in  megahertz) 

receiver  noise  bandwidth  (in  megahertz) 

receiver  noise  figure  (in  decibels) 

integration  improvement  in  S/N  ratio  (in  decibels)  (e.g. , 
if  the  radar  integrates  100  pulses,  perfect  integration 
improvement  would  be  entered  as  20)  (see  reference  1, 
pages  35-40) 

radar  system  losses  (in  decibels) 

radar  pulse  length  (in  microseconds) 

azimuth  beamwidth  in  degrees  measured  at  the  3 db 
down  level 

radar  antenna  hei^t  (in  feet) 

radar  polarization  (0=vertical;  l=horizontal) 

output  in  graph  form  (l=yes;  0=no) 


Effective  Earth's  Radius  « X(15)  times  earth's  radius 


Using  X(15)  = 1.333  corresponds  to  the  nominal  4/3  earth  approximation 


Atmospheric  attenuation,  the  absorption  of  a portion  of  the  radar's  signal  strength  by 
the  atmosphere,  will  be  discussed  in  a later  section.  The  technique  this  detection 
model  uses  to  calculate  the  amount  of  attenuation  was  determined  from  attenuation  curves 
contained  in  reference  5.  These  attenuation  curves  assume  the  Central  Radio  Propagation 
Laboratory  exponential  reference  atmosphere  for  refraction  and  the  International  Civil 
Aviation  Organization  standard  atmosphere  for  p'^.isure-temperature  values.  The  detec- 
tion model  cannot  consider  irreguiar  atmospheres. 


Noise  jamming  signals  are  enhanced  by  the  radar's  antenna  gain  on  the  Jamming 
bearing.  Figure  2 presents  an  example  of  a radar  antenna  gain  pattern  with  the  main-lobe 
and  various  side-lobes  described  by  the  shaded  area.  This  model  uses  a single  value 
(averaged  throughout  the  side-lobes)  for  antenna  gain  in  the  side-lobes.  This  method  is 
represented  in  figure  2 by  the  dashed  line.  A closer  approximation  to  reality  could  be 
achieved  by  allowing  the  jammers  to  operate  against  the  antenna  gain  averaged  over  a 
smaller  region.  This  technique  is  represented  by  the  dotted  line  in  figure  2. 


FIG.  2:  SAMPLE  ANTENNA  GAIN  PATTERN 
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The  model  requires  that  all  the  staodoff  jammers,  both  In  the  main-  and  side -lobes, 
be  at  the  same  range  and  have  identical  jamming  bandwldtha. 

The  model  does  not  consider  the  multipath  effects  of  the  jamming  signal  Interfering 
with  Itself,  nor  does  it  allow  for  the  movement  of  the  jamming  vehicle. 

CAM  DETECTION  MODE  I 

The  CAM  detection  model  is  the  basis  for  this  revised  model.  The  model  calculates 
the  signal-to-nolse  ratio  as  a function  of  target  range  using  the  following  equations: 


^^^TAR  ^MUL2  ^ATT2 

^REC  ^ATT2  ^^^‘^SCLTR  ^^^WTHR^"^  ^ATTl^^^JAM 


where: 


Fatti  “ one-way  attenuation  factor 
Fatt2  ~ two-way  attenuation  factor 
o “ two-way  multipath  effect 
= signal  received  from  the  target 


Nrec  = receiver  noise 

SIGgcLTR  “ signal  received  from  surface  clutter 

SIC™—,!,  = signal  received  from  weather  clutter 
W 1 HK 

SIGjam  = signal  received  from  the  jamming  source 
The  signal  received  from  the  target  is  calculated  using  equation  2: 


(4lT)3LgR^< 


where: 


= peak  transmitted  power  (watts) 

G.J,  = antenna  gain,  main-lobe  (decibels) 

ne^  (n)  = integration  improvement  in  S/N  ratio  (decibels) 


X = wave  length  (meters) 

= average  target  cross  section  (square  meters) 
= range  from  radar  to  target  (meters) 

L-  = system  losses  (decibels) 

d 

The  receiver  noise  is  as  in  equation  3: 


^REC  ^ 


where: 

-23 

K = Boltzmann's  constant  (1.38X10  Joules /degree) 

Tq  = standard  temperature  (290°  K) 

= receiver  noise  figure  (decibels) 

Bv,  = receiver  noise  bandwidth  (hertz) 

N 

The  signal  received  from  surface  clutter  is  calculated  using  equation  4: 


SCLTR 


^tVx^B  CTO 
Z o 

3 3 
2(4tt)  R 


where: 


P.J.  = peak  transmitted  power  (watts) 

G^=  antenna  gain,  main-lobe  (decibels) 

X = wave  length  (meters) 

B„  = azimuth  bandwidth  of  radar  (degrees) 

" 8 
C = velocity  of  light  (2.997925  x 10  meters  per  second) 

T = pulse  length  (seconds) 

o^  = surface  clutter  coefficient  (decibels) 

R = range  at  which  S/N  ratio  is  being  calculated  (meters) 


The  signal  received  from  weather  clutter  Is  written  as  follows: 

.93  P_G_CtttS 

“ 128  X2  <5) 

where: 

= peak  transmitted  power  (watts) 

G— = antenna  gain,  main*lobe  (decibels) 

^ g 

C = velocity  of  light  (2.997925  X 10  meters  per  second) 

T = pulse  length  (seconds) 

X = wave  length  (meters) 

R = range  at  which  S/N  ratio  is  being  calculated  (meters) 

Z = 200  r^*^  (10  ^®)  (cubic  meters) 
r = rainfall  rate  (millimeters  per  hour) 

The  signal  received  from  the  jammer  is  calculated  according  to  equation  6: 


8IG 


JAM 


2 

(4„)'= 


(6) 


where: 


X 

B. 


effective  radiated  power  of  noise  jammer  (watts) 

effective  radar  gain  in  the  direction  of  the  jammer  (main-  or  side-lobe) 
(decibels) 

wave  length  (meters) 
receiver  noise  bandwidth  (hertz) 
noise  jammer  bandwidth  (hertz) 
range  from  radar  to  jammer  (meters) 


The  signal  received  from  the  jammer  must  be  modified  by  a one-way  atmospheric 
attenuation  factor: 


ATT 


(7) 


-9- 


where: 


a = attenuation  coefficient  (meters  S ^ 

Rj=  range  from  radar  to  jammer  (meters) 

The  signals  received  from  the  target,  surface  clutter,  and  weather  clutter  must  be 
modified  by  a two-way  atmospheric  attenuation  Victor: 

„ -2a  R 

^ATT2~  ® 

where: 

a = attenuation  coefficient  (meters 

R=  range  at  which  S/N  ratio  is  being  calculated  (meters) 

Tne  signal  received  from  the  target  must  be  modified  by  the  two-way  multipath  propa- 
gation factor  (reference  4): 

'’muli  ' ^ 

F = (F 
MUL2  ' MULl^ 

where: 

^MULl  ~ multipath  propagation  factor 

F 

MUL2  = two-way  multipath  propagation  factor 

P = power  at  the  target  including  multipath 

Pq  = power  at  the  target  excluding  multipath 

R^  = smooth  sea  reflection  coefficient 

p = rough  sea  reflection  coefficient 

D = divergence  factor 

= phase  angle  of  direct  ray  relative  to  reflected  ray 

There  are  two  common  uses  of  the  symbol  "o-"  within  the  radar  community.  The  first 
usage  refers  to  the  1 W/MHZ/M^  self -screening  jammer  detection  range  and  the  second 
refers  to  the  atmospheric  attenuation  coefficient. 

-10- 


(9a) 

(9b) 


The  smooth  sea  reflection  coefflecient  of  the  reflected  ray  relative  to  the  incident  ray  is 
given  by  (reference  4): 


(vertically  polarized  radar) 


(horizontally  polarized  radar)  (11) 


Ry  = magnitude  of  the  reflected  ray  (vertical  polarization) 

= magnitude  of  the  reflected  ray  (horizontal  polarization) 

= phase  angle  of  the  reflected  ray  (vertical  polarization) 

= phase  angle  of  the  reflected  ray  (horizontal  polarization) 
n = complex  dielectric  constant 
X = wave  length 

= permitivlty  constant  of  the  reflecting  surface 
a = conductivity  constant  of  the  reflecting  sur&ce 
Y = angle  of  Incidence  = angle  of  reflection 


For  rough  seas,  the  magnitude  of  the  reflected  ray  relative  to  the  incident  ray  is  given  by 
(references  4 and  5): 


where; 


h = root -mean -square  wave  height 
Y = angle  of  incidence  = angle  of  reflection 
X = wave  length 


The  phase  shift  caused  by  path  length  difference  is  given  by  (reference  2): 


2n 

X 


(13) 


where: 

AR  = path  length  difference 
X = wave  leitgth 


The  total  phase  angle  is  given  by  (reference  4): 

(vertical  polarization) 

4>  = ^5  - or 

h (horizontal  polarization) 


(14) 


The  divergence  factor  is  a result  of  the  ray  being  reflected  from  a spherical  surface 
rather  than  a flat  jurface.  Divergence  spreads  the  ray,  reducing  the  signal  strength 
per  unit  area;  and  thereby  reducing  the  multipath  effect.  The  divergence  factor  is  given 
by  (reference  4): 


D 


R (R,+R-)sin 
e X ^ 


i-n  yJ 


(15) 


SELF -SCREENING  JAMMER 

The  Q and  0 ranges  are  commonly  used  measures  of  radar  performance; 

Q range  — The  detection  range  against  a one  square  meter  target  with  a one 
watt  per  megahertz  self -screening  jammer. 

jS  range  — The  detection  range  against  a one  square  meter  target  in  a clear 
ECM  environment. 

A radar  detection  model  can  demonstrate  to  the  analyst  that  the  radar  is  being  mode.'ed 
correctly  by  comparing  the  a and  0 ranges  to  values  obtained  from  other  reliable 
sources . 

By  use  of  input  ”X(10), " a self -screening  jammer  capability  is  treated.  This  input 
is  programmed  such  that,  if: 

X(10>=0  The  model  would  accept  a standoff  jammer  with  the  jammer's 
range  being  the  value  input  as  X(12) 

X(10)=l  The  model  would  accept  a self -screening  jammer  with  the 
jammer's  range  being  the  target's  range 


-12- 
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ASCENDING/DESCENDING  TARGETS 

The  target's  position  at  tiie  beginning  and  end  of  the  run  is  entered  in  relation  to  the 
surface  of  the  earth  (see  figure  3).  These  positions  are  converted  into  polar  coordinates. 

^ _ TT  RNGBGN 
1~2  " RE+ALTBGN 

where; 

Gj  = target's  angular  (radians)  position  at  the  beginning  of  the  run 
RNGBGN  = target's  range  (feet)  at  the  beginning  of  the  run 
RE  = radius  of  the  earth  (feet) 

ALTBGN  = target's  altitude  (feet)  at  the  beginning  of  the  run 

The  target's  position  in  polar  coordinates  at  the  beginning  is  (PREALT,  G^)  where 
present  altitude  (PREALT)  is  equal  to  the  sum  of  RE  and  ALTBGN. 


The  target’s  position  at  the  end  of  the  run  is  determined  by: 

^ ^ rr  RNGEND 
2 2 " RE  + ALTEND 

where: 

^2  = target’s  angular  position  (radians)  at  the  end  of  the  run 
RNGEND  = target’s  range  (feet)  at  the  end  of  the  run 
ATTEND  = target’s  altitude  (feet)  at  the  end  of  the  run 

The  target’s  end  position  in  polar  coordinates  is  (FNLALT,  D ) where  the  final  altitude 
(FNLALT)  is  equal  to  the  sum  of  RE  and  ATTEND. 

The  target’s  rate  of  ascent/descent  (y)  is  determined  by: 

ALTBGN-ALTEND 

CWT • 08) 

2 1 


The  model  calculates  the  S/N  ratio  at  incremental  ranges  from  the  beginning  through 
the  end  of  the  run.  The  incremental  change  is  determined  fay: 


STEP  = 


1 


RNGBGN-RNGEND 


(19) 


During  the  first  iteration,  the  model  converts  the  beginning  position  (PREALT,  0^) 
into  rectangular  coordinates,  and  then  calculates  the  S/N  ratio. 


HT=  PREALT  - RE  (target  height  in  feet) 


RT  = '^[(sin  flj)  (PREALT)  - (RE+H^)]  ^ + [(cos  ^2)  (PREALT)]^  (20) 

where: 

RT  = target  range  (feet) 

= radar  antenna  height  (feet) 
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Following  this  calculation,  the  model  decreases  the  target  range  by  one  incremental 
step. 


-*■  STEP  (new  angular  position  of  target)  (21) 

PREALT  + PREALT  - [(y)  (STEP)]  (new  altitude  of  target)  (22) 


The  model  again  converts  the  target’s  position  into  rectangular  coordinates  and 
calculates  the  S/N  ratio.  The  model  continues  to  close  the  target  toward  the  radar  in 
incremental  stq;)s  until  the  target  reaches  its  end  position. 

If  the  target  has  more  than  one  rate  of  ascent  or  descent  during  its  flight  profile, 
the  rates  must  be  separated  into  individual  cases  involving  single  rates  of  ascent  or 
descent.  The  model  can  also  allow  for  targets  at  a constant  altitude  by  entering  both 
the  beginning  and  ending  altitude  as  the  same  value. 

SIMULTANEOUS  MAIN-LOBE  AND  SIDE -LOBE  JAMMING 

With  a multi -axis  threat,  one  might  expect  simultaneous  main-lobe  and  side-lobe 
jamming. 


The  location  of  the  jamming  is  programmed  using  input  X(ll),  if: 


X(ll)  = 0 side -lobe  jamming  only 

X(ll)  = 1 main-lobe  jamming  only 

X(ll)  = 2 both  main-lobe  and  side-lobe  jamming. 

The  power  of  the  jammer(s)  is  programmed  using  input  X(14). 

X(14)  = SSSSMMMM, 

a packed  word  with  "SSSS"  referring  to  side -lobe  jamming  power  and  "MMMM"  being 
main-lobe  jamming  power,  in  watts. 


The  signal  received  from  the  jammers  (see  equation  6)  is  separated  into  two  parts: 


SIG 


JAM(M) 


■ ‘’kM)  °KM)  °'nM)  * ^ 
(41^*  B|  Rj  * 


(main-lobe  jamming)  (23) 
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SIG 


JAM  (S) 


- ’’kS)  °1(M)  °TtS)  ^ °N 


2 2 
(4tt)^  Rj 


(side -lobe  jamming)  (24) 


where: 

^J(M)  ‘^KM) 
°T(M) 

^](S)  ^J(S) 
°T(S) 


effective  power  of  all  Jammers  radiating  in  the  main -lobe 
effective  main-lobe  radar  gain 

effective  power  of  all  jammers  radiating  in  the  side -lobe 
effective  side -lobe  radar  gain. 


ATTENUATION 

The  atmospheric  attenuation  of  a radar's  signal  strength  in  a clear  atmosphere  is 
due  primarily  to  oxygen  and  water  vapor.  A portion  of  the  radiated  energy  is  absorbed 
as  heat  by  these  atmospheric  factors  and  is  lost. 


Figures  4 and  5 present  attenuation  curves  taken  from  reference  5.  These  curves 
are  treated  as  a function  of  target  range,  radar  frequency,  and  elevation  (horizontal  to 
target)  angle.  These  curves  were  calculated  using  the  Central  Radio  Propagation  Labora- 
tory exponential  reference  atmosphere  for  refraction  and  the  International  Civil  Aviation 
Organization  standard  atmosphere  for  pressure -temperature  values.  This  section 
describes  the  technique  to  calculate  attenuation,  which  is  based  on  the  attenuation  curves 
from  reference  5.  Irregular  atmospheres  cannot  be  considered  in  this  detection  model. 


The  CAM  detection  model,  as  well  as  references  1 and  5,  treats  the  effects  of 
atmospheric  attenuation  by  an  exponential  law. 


F 

ATTl 


(one-way  attenuation  factor) 


F 

ATT2 


-2oR 
= e 


(two-way  attenuation  factor) 


where  a is  the  attenuation  coefficient  and  R is  the  target  range.  The  reader  should  note 
that  the  attenuation  coefficient  is  multiplied  by  two  to  obtain  two-way  attenuation.  The 
attenuation  coefficient  is  determined  by  dividing  the  one-way  attenuation  loss  by  the 
target  range. 

The  attenuation  coefficients  for  a 1000  megahertz  radar  with  a OP  elevation  angle 
were  calculated  at  20  nautical  mile  intervals  from  figure  4 and  graphed  on  figure  6.  The 
curve  of  the  attenuation  coefficients  is  not  linear  with  range.  Using  a constant  value  for 
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Radar  frequency  (Mhz)  10,00 


FIG.  4:  ABSORPTION  LOSS  FOR  TWO-WAY  PROPAGATION  AS  A FUNCTION  OF 
RADAR  RANGE,  VARIOUS  FREOUENCIES,  AND  OP  ELEVATION  ANGLE 
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FIG.  5:  ABSORPTION  LOSS  FOR  TWO-WAY  PROPAGATION  AS  A FUNCTION  OF  RADAR 
RANGE,  VARIOUS  FREQUENCIES,  AND  0.5°  ELEVATION  ANGLE 


COEFFICIENT  OF  ATTENUATION  FOR  A RADAR  WITH  A FREQUENCY  OF 
1.000  MEGAHERTZ  AND  A 0°  ELEVATION  ANGLE 
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the  attenuation  coefficient  introduces  an  error  that  can  be  significant  for  high  frequency 
radar  at  low  elevation  angles. 

The  following  technique  was  devised  to  reduce  the  errors . Users  of  this  program 
will  need  to  use  this  technique  for  the  particular  frequency  and  elevation  angle  of  the 
radar  they  plan  to  examine.  For  a 2-D  radar,  only  one  calculation  nuiy  need  to  be 
made  at  the  elevation  angle  of  the  radar's  main  beam.  For  3-D  radars,  calculations 
should  be  made  for  each  beam  grouping.  Again,  for  3-D  radars,  the  program  should 
be  run  as  separate  cases  corresponding  to  the  target's  flight  profile  through  the  appro- 
priate beam  grouping  while  using  the  attenuation  for  that  beam  grouping. 

Table  2 presents  the  one-way  attenuation  losses  for  a 1000  megahertz  radar  and 
0°  elevation  angle. 

TABLE  2 


ATTENUATION  LOSS 


Range 

(n.mi.) 

20 

40 

60 

80 

100 

120  140  160 

180  200 

220  240 

260  280  300 

Loss 

(db) 

.21 

.40 

.59 

.77 

.93 

1.08  1.20  1.29 

1.35  1.40 

1.45  1.48 

1.49  1.50  1.50 

The  values  were  determined  by  taking  the  appropriate  values  from  figure  4 (two-way 
attenuation)  and  dividing  them  by  two.  A natural  logarithmic  curve-fitting  technique 
is  used  to  determine  an  attenuation  coefficient.  The  following  regression  equation  is 
used: 

or  = Oj  + R (25) 

To  simplify  the  work,  the  regression  coefficients,  and  , are  changed  to  attenua- 
tion losses  from  the  atmospheric  attenuation  curve  (figure  4)  by  first  multiplying  equation 
25  by  range: 


Ra  = Ro,  + Rof.  Ln  R 

1 2 

(26a) 

L = L,  + L„  Ln  R 

1 2 

(26b) 
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where  L is  the  attenuation  loss  corresponding  to  range  R (these  values  are  identified 
in  table  2).  A regression  is  accomplished  and  and  are  inserted  into  the  radar 

detection  model  as  inputs  X(7)  and  X(8).  The  radar  detection  model  will  take  the  inputs, 
calculate  the  attenuation  ioss  at  the  target  range  being  examined,  convert  the  attenuation 
loss  to  an  attenuation  coefficient  in  decibels  per  nautical  mile,  and  convert  it  to  an 
attenuation  coefficient  per  nautical  mile. 

As  an  example,  a natural  logarithmic  regression  on  the  data  in  table  2 was  accom- 
plished with  the  ioss  coefficients  determined  to  be: 

L = Ln  R (see  equation  26b) 

L=  - 1.55605  + 0.54851  Ln  R . 


A correlation  coefficient  of  0.98  was  determined  for  the  regression.  Using  the  t-test, 
the  correlation  coefficient  tested  significant  at  the  99  percent  confidence  level. 

The  results  of  using  the  natural  logarithmic  regression  method  are  displayed  in 
figure  7.  The  theoretical  loss  curve  was  taken  from  figure  4. 

In  a later  section,  a fictitious  radar  will  be  modeled  for  example  purposes.  The 
regression  coefficients  determined  above  will  be  used  as  inputs  to  the  detection  modei. 

MULTIPATH 

The  effect  of  multipath  interference  occurs  when  the  radar  signal  reflects  off  the 
surface  such  that  some  of  the  radiation  reaches  the  target  by  a reflected  path  (R^  and  R^, 

see  figure  5)  in  addition  to  the  direct  path  (RT).  These  two  rays  generally  arrive  at  the 
target  with  both  a phase  and  amplitude  difference.  The  field  strength  is  the  phasor  sum 
of  these  two  rays. 

The  geometry  of  the  multipath  problem  is  shown  in  figure  8,  where: 

RT  = target  range 

= radar  antenna  height 
= target  altitude 

RE  = effective  radius  of  the  earth  that  includes  refraction 
^^=^2  = angle  of  incidence  = angle  of  reflection 
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FIG. 


Reference  4,  using  the  assumptions: 


RT  *»  long  range  with  Y small 


-;-2RE — >0 


shows  that: 


T,  3 

r 3 o_i 

2 fl  2 1 

r 1 

Rl  + 

[-T-J 

Rl  + 2 -RE(Hj+H2)  Rj+ 

RE  • Hj  • RT 
(29) 

Reference  4 then  presents  a method  to  determine  which  one  of  the  three  roots  produces  the 


actual  value  for  R^  . 

Following  the  calculation  of  Rj,  reference  4 uses: 


R2  = H2  + sin  Y 


to  determine  the  value  of  R^.  The  accuracy  of  R2  decreases  at  long  ranges  using  this 

equation.  From  the  geometry  of  figure  8 it  can  be  seen  that  the  side  of  the  triangle  is 
not  H2  but  minus  X^. 

Computation  Method 

Although  not  mathematically  pleasing,  the  following  method,  which  has  been  inserted 
into  the  revised  radar  detection  model,  will  calculate  the  path  length  difference  with  an 
error  of  less  than  one  inch  out  to  a target  range  of  300  nautical  miles.  It  takes  31  iterations 
to  produce  such  a level  of  accuracy.  The  method  used  requires  no  assumptions  and  is 
therefore  useful  at  all  target  ranges  and  altitudes . 

From  the  Law  of  Cosines  it  can  be  shown  that: 


cos  (9j+02)  = 


(RE+H^)^-KRE+H2)^-RT^ 


2 (RE+Hp  (RE+H2) 


0^+02  = arc  cos 


cos  (9i+02) 


Although  the  sum  of  6^  plus  ^ found,  the  individual  values  are  not  known. 

Define  an  interval  0 degrees  to  degrees.  In  the  first  iteration,  let  0j=l/2  (®2+02)* 

Using  this  assumption,  can  be  computed  from  the  Law  of  Cosines: 


RE^  + (RE+H  )^-R  ^ 
® 1 2 • RE  • (RE+H,) 


(33) 


Equation  33  can  be  solved  for  R^: 


RE"  + (RE+Hj)  -2(RE)  (RE+H^)  cos  0^ 


(34) 


! S 


Using  the  R,  determined  by  the  assumed  0j,  can  be  computed  using  the  Law  of 
Cosines: 


Similarly, 


^ RE^+R  ^-(RE+H  )^ 

cos  (j+  Yj)  = 2(RE)(R,) 


-sin  = cos  (j+  'i’j) 


(trigonometric  identity) 


sin  = 


(RE+Hj)2  -(RE^)  -(R^2) 
2(RE)(Rp 


= arc  sin  (sin  Y^) 


R = \f\ 


RE^-+<RE+H2)^-2(RE)  (RE+H^)  cos  02 


^2  = arc  sin 


(RE+H2)^-RE^  -R2^ 


2(RE)(R2) 


(35) 


(36) 


(37) 


(38) 


Now  test  against  . If  'l'j>'f2  » redefine  the  interval  as  0^  degrees  to  ( ^2+^2) 
degrees.  If  Yj<Y2  , redefine  the  interval  as  0 degrees  to  0^  degrees.  Recalculate 
Yj^  and  with  0j  assumed  to  be  one-half  of  the  new  Interval. 
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The  program  in  appendix  E was  taken  from  the  revised  detection  model  to  test  the 
accuracy  of  the  method.  Tables  E-1  through  E-4  are  results  at  various  target  ranges 
and  altitudes . 

An  error  of  les^  thanj)ne  inch  in  the  point  of  reflection  would  create  a negligible 
error  in  the  phase  shift  caused  by  the  path  length  difference. 


Oj 


This  method  will  reduce  the  error  between  the  initially  assumed  0 and  the  actual 

N ^ 

0j  by  a factor  of  1/2  (where  N is  the  number  of  iterations),  fa  31  iterations,  this 

method  will  reduce  the  error  between  and  ®“ch  that  the  error  between  the  cal- 
culated point  of  reflection  and  the  actual  point  of  reflection  is  less  than  one  inch. 


Oj 


Oi 


^ 1 
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APPENDIX  A 
USERS’  GUIDE 


i t 


This  appendix  provides  users*  instructions  for  the  radar  detection  model.  The 
control  cards,  the  DATA2  subroutine,  and  the  PLOTTER2  subroutine  are  specifically 
for  CNA’s  computer  (GDC  3800).  It  is  reconunended  that  non-CNA  users  replace  the 
DATA2  and  PLOTTER2  subroutines  with  similar  routines  within  your  own  system.  The 
radar  detection  model  is  written  in  FORTRAN,  which  may  be  useable  on  other  com- 
puters with  minimal  modifications. 


Control  Cards 


A tjpical  card  deck  is  shown  in  figure  A-1  . The  DATA2  subroutine  for  reading 
in  data  allows  for  several  cases  to  be  run  at  one  time. 


DATA2  Subroutine 


The  DATA2  subroutine  is  used  to  read  in  34  input  parameters  and  locate  them  in 
an  array.  The  subroutine  has  two  beneficial  features.  It  allows  free  field,  formatless 
use  of  input  cards,  and  the  contents  of  the  array  are  zeroed  only  at  the  beginning  of  a 
job,  not  between  cases. 
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In  using  the  free  field  feature,  the  first  number  on  the  input  card  identifies  an  ele- 
ment in  the  input  array.  The  next  number  is  the  value  that  will  be  placed  into  that  I 

element  specified  by  the  first  number.  The  third  number  is  the  value  that  will  be  placed  ^ 

in  the  next  sequential  element  of  the  array.  Additional  numbers  will  be  placed  in  sub-  I 

sequent,  sequential  elements.  There  must  be  at  least  one  space  between  each  number  j 

(e.g.,  78  135  1E6. . .).  Input  cards  are  read  in  sequence  for  a given  case  until  a blank  I 

card  is  reached.  ^ 

i 

The  array  is  zeroed  at  the  beginning  of  a job.  All  values  used  in  the  first  case  of  I 

a job  are  continued  into  the  second  case,  except  those  specifically  changed.  The  same  i 

process  occurs  for  subsequent  cases.  This  feature  allows  changing  parameters  between 
cases  without  having  to  reenter  the  entire  input  deck.  Only  those  parameters  to  be 
changed  need  be  specified  after  the  first  case.  Documentation  and  a more  thorough 
discussion  of  DATA2  are  found  in  reference  A -6.  i 

The  following  example  may  assist  the  user  in  understanding  the  use  of  DATA2.  i 

Figure  A-?,presents  input  decks.  Sample  case  1 has  six  values  (table  A-1)  to  be  i 

entered.  ' j 


FIG.  A-1:  TYPICAL  CARD  DECK 
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TABLE  A-1 

f 

INPUT  VALUES 

Element  in  the  array 

Vail 

1 

16 

if 

2 

80 

3 

.5 

4 

16 

5 

-12 

f .. 

6 

.0 

FIG.  A-2:  DATA2  EXAMPLE 


In  sample  case  2,  we  wish  to  change  input  4 from  16,  000  to  17,  000.  Note  that  the  other 
values  need  not  be  included.  In  sample  case  3,  we  wish  to  change  input  5 from  -12  to  -10 
and  input  6 from  .001  to  .002  . Again,  inputs  that  have  not  been  changed  remain  the  same 
as  in  the  previous  case.  Table  A-2  Identifies  the  input  values  to  be  used  in  sample  case  3. 
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TABLE  A-2 
INPUT  VALUES 


Value 

16 

80 

.5 

17,  000 
-10 
.002 


PLOTTER2  Subroutine 

The  PLOTTER2  subroutine  takes  the  output  from  the  radar  detection  model  and 
plots  the  signal-to-noise  ratio  as  a function  of  target  range.  Documentation  of  PLOTTER2 
is  contained  in  reference  A-7. 


If  the  program  is  to  be  run  on  CNA’s  computer,  the  user  need  not  be  concerned  with 
this  subroutine.  Programming  has  been  inserted  to  make  the  plotting  automatic.  The 
plot  will  be  7 x 10  inches  (vertical  x horizontal).  TTie  output  is  scaled  for  maximum 
projection  within  the  7 x 10  inch  limits.  For  example,  the  horizontal  axis  has  10 
scaling  marks.  If  the  maximum  target  range  is  50  nautical  miles,  the  scaling  marks 
will  represent  5 nautical  miles.  If  the  maximum  target  range  is  150  nautical  miles, 
the  scaling  marks  will  represent  15  nautical  miles.  See  appendixes  Band  C for  ex- 
amples of  the  output  plots. 
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APPENDIX  B 


STANDOFF  JAMMING  EXAMPLE 

In  order  to  assist  a future  user  to  understand  the  radar  detection  model,  two 
examples  will  be  provided.  It  is  emphasized  that  the  radar  and  the  threats  in  the 
scenario  are  fictitious.  The  two  examples  will  use  the  same  radar,  threat  missile, 
and  environmental  conditions.  In  this  appendix,  we  shall  examine  standoff  noise 
jamming.  The  following  appendix  will  assume  a missile  with  a self-screening  noise 
jammer  but  no  standoff  jamming. 

The  basic  scenario  is  as  follows.  The  ship  is  in  sea  state  two  (Beaufort  scale) 
and  multipath  interference  can  be  expected.  The  atmospheric  attenuation  can  be  des- 
cribed by  the  regression  coefficients  calculated  in  the  previous  section  on  attenuation. 

L — Lj^  + Ij^  Ln  R 

L - 1.55605  +0.54851  Ln  R 

A surface  clutter  coefficient  of  -30  decibels  was  determined  from  reference  1 (pages 
527-534).  There  is  a drizzle  with  a rainfall  rate  of  1 millimeter  per  hour.  A coeffic- 
ient of  refractivity  of  1.333  is  assumed,  corresponding  to  the  4/3  earth  approximation. 

The  user  has  researched  various  reference  publications.  Table  B-I  lists  the 
characteristics  of  the  radar  determined  from  the  references. 

The  fictitious  threat  missile  has  a one  square  meter  radar  cross  section.  The 
missile  is  launched  from  150  nautical  miles  at  an  altitude  of  15, 000  feet.  The  missile 
has  an  unspecified  homing  device  that  is  locked  onto  the  ship  prior  to  launch.  The 
missile  flies  a constant  azimuth,  constantly  decreasing  altitude  flight  profile  from  launch 
at  15, 000  feet  until  it  hits  the  ship. 

There  are  four  enemy  jamming  aircraft  at  150  nautical  miles.  Each  jamming  air- 
craft is  separated  by  90  degrees.  Each  aircraft  has  a jamming  capability  of  1,  000  watts 
over  a bandwidth  of  10  megahertz.  The  threat  missile  may  be  launched  from  any  one 
of  the  four  aircraft.  As  a result,  the  threat  missile  will  always  be  in  a main-lobe 
jamming  strobe  of  one  of  the  four  aircraft.  The  other  three  aircraft  will  be  jamming 
into  the  radar's  side-lobes. 

The  inputs  to  the  detection  model  are  presented  in  table  B-2.  Figure  B-1  provides 
the  input  deck.  The  numerical  outputs  (table  B-3)  and  the  graphical  outputs  (figure  B-2) 
present  the  results  of  the  first  example. 


B-1 


TABLE  B-1 


FICTITIOUS  RADAR  PARAMETERS 


Parameter 

Peak  transmitted  power  (kw) 

Antenna  gain,  main-lobe  (db) 

Antenna  gain,  side-lobe  (db  down  from 
main- lobe,  db) 

Operating  frequency  (Mhz) 

Receiver  noise  bandwidth  (Mhz) 

Receiver  noise  figure  (db) 

Integration  improvement  in  S/N  ratio  (db) 

System  losses  (db) 

Pulse  length  (usee) 

Azimuth  beamwidth  (down  3 db  level,  deg? 
Antenna  height  (ft) 

Polarization 


16  2 1SOOO  0 150  0 1 1000  35  25  1000  10  5 50  1 0.5  5 100  1 1 


1 1 1 -30  1 1 1 -1.55605  0.54851  1 0 2 150  10  30001000  1.333 


Value 


1000 


1000 


Horizontal 


FIG.  B-1:  STANDOFF  JAMMING  EXAMPLE  INPUT  DECK 
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TABLE  B-2 

RADAR  DETECTION  MODEL  INPUTS 


IHPOT 


DBFINITION 


VALUE 


■ultipath  effect  included  (l>yes;  0«no) 

surface  clutter  Included  (l-yes;  0-no) 

surface  clutter  coefficient  in  decibels. 

A ratio  of  1 square  meter  of  clutter  to 
1,000  square  meters  illuminated  would  be 
entered  as  -30  (see  reference  1,  pages 
527-534) 

weather  clutter  Included  (1-yes;  0-no) 

rainfall  rate  in  millimeters  per  hour 

attenuation  Included  (1-yes;  0-no) 

X(7)-L]^  and  X(8)-L2  of  natural  logarithmic 
regression  L-L^^+L^LnR,  where  L is  the 
attenuation  loss  in  decibels  and  R is 
target  range  in  nautical  miles  (see  section 
on  attenuation 


0.54851 


jamming  included  (1-yes;  0-no) 

type  of  jammer  (0-standoff  jammer; 
1-self-screening  jammer) 


X(10) 


X(ll) 


location  of  jammming  (0-side-lobe; 

1-main- lobe;  2-both  main-  and  side-lobes) 


X(12) 


range  from  radar  to  jammers  in  nautical  ^50 

miles.  Note:  all  j2unmers  must  be  at 
the  same  range 

noise  jammer  bandwidth  in  megahertz.  Note: 
all  j2unmers  must  have  identical  bandwidths 

effective  radiated  power  of  the  jammer(s)  03000100 
in  watts.  Pa^;ked  word  "SSSSMMHM"  where 
•SSSS"  refers  to  side-lobe  jamming  and 

refers  to  main-lobe  jamming.  Note: 
all  jamming  power  in  the  appropriate  lol>e 
must  be  totaled  and  entered  as  one  value. 

For  example,  three  100  watt  jammers  in  the 
side-lobe  and  two  100  watt  jammers  in  the  main- 
lobe  would  be  entered  as  03000200  . 


X(13) 


X(14) 


Ti\BLE  B-2  (Cont’d) 


INTOT 

DEFINITION 

VALUE 

X(15) 

coefficient  of  refractivlty  (1.333  for  the 

4/3  earth  approxiaatlon) 

1.333 

X(16) 

aea  state  (Beaufort  scale) 

2 

X(17) 

target  altitude  at  the  beginning  of  the 
run  (In  feet) 

15,000 

X(18)  . 

target  altitude  at  the  end  of  the  run  (In 
feet) 

0 

X(19) 

target  range  at  the  beginning  of  the  run 
(in  nautical  miles) 

150 

X(20) 

target  range  at  the  end  of  the  run  (in 
nautical  miles) 

0 

X(21) 

average  target  cross  section  (in  square 
SMters) 

1 

X(22) 

Pealc  radar  transmitted  power  (in  Icilowatts) 

1000 

X(23) 

main-lobe  antenna  gain  (in  decibels) 

35 

X(24) 

side-lobe  antenna  gain  (in  decibels) (e.g. , 
side-lobe  down  2S  db  from  the  main-lobe 
gain  would  be  entered  as  25) 

25 

X(2S> 

radar  frequency  (in  megahertz) 

1000 

X(26) 

receiver  noise  bandwidth  (in  megahertz) 

10 

X(27) 

receiver  noise  figure  (in  decibels) 

5 

X(28) 

integration  improvement  in  S/N  ratio  (in 
decibels) (e. g. , if  the  radar  integrates  100 
pulses,  perfect  integration  improvement  would 
be  entered  as  20) (see  reference  i.  pages 
35-40) 

50 

X(29) 

radar  system  losses  (in  decibels) 

1 

X(30) 

radar  pulse  length  (in  microseconds) 

.5 

X(31) 

azimuth  beamwidth  in  degrees  measured  at 
the  3 db  down  level 

5 

X(32) 

radar  antenna  height  (in  feet) 

100 

X(33) 

radar  polarization  (0-vertical;  1-horizontal) 

1 

X(34) 

output  in  graph  form  (1-yes;  0-no) 

1 
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TABLE  B-3 

STANDOFF  JAMMING  EXAMPLE 


0.0 

05.2 

50.0 

32.8 

100.0 

8.3 

1.0 

"127.6 

51.0 

31.0 

101.0 

9.4 

3.0 

2.7 

52.0 

22.0 

102.0 

6.5 

3.0 

36.3 

53.0 

27.1 

103.0 

4.6 

M.O 

40.2 

54.0 

22.0 

104.0 

1.8 

5.0 

54.8 

55.0 

10.1 

105.0 

"11.0 

6.0 

57.8 

56.0 

22.6 

106.0 

"7.3 

7.0 

50.1 

57.0 

17.3 

107.0 

"5.0 

8.0 

50.6 

58.0 

5.7 

108.0 

"19.0 

0.0 

18.2 

50.0 

15.7 

109.0 

"19.1 

10.0 

55.1 

60.0 

"10.7 

110.0 

"11.4 

11.0 

44.3  • 

61.0 

"28.7 

111.0 

"22.4 

12.0 

40.5 

62.0 

"6.0 

112.0 

"14.3 

13.0 

24.1 

63.0 

_^"4.5 

113.0 

"1.1 

14.0 

15.0 

7.5 
31.9  / 

-N 

64.0 

65.0 

"38.2 
/ 4.5 

1T470">' 

~115.0/ 

"1.1 

"0.0 

16.0 

17.0 

40.2/ 

4.0 

ii 

N \ 

Ni 

66.0 

67.0 

0.8  - 
, 10.6 

116-.-0 

llT^.O 

1.5 

7.5 

18.0 

21.6 

60}o 

/ 10.6 

ilO.O 

7.6 

10.0 

27.3, 

\ '■ 

60  ;o 

18.2 

. '119.0 

8.8 

20.0 

"27.2  V 

' V. 

70.0 

13.7 

X 120.0 

11.7 

21.0 

26.7 

' - 

71.0 

'16.1 

' 121.0 

12.1 

22.0 

"5.5 

/ , 

72.0 

21.8 

/ 122.0 

12.5 

23.0 

27.4 

/ 

73.0 

21 . 5 

123.0 

14.1 

24.0 

"32.1 

/ 

1 

74.0 

22.1' 

124.0 

13.8 

25.0 

26.5 

✓ 

/ 

/ 

75. 0 

24.2 

125.0 

14.4 

26.0 

12.0 

\76.0 

.21.3" 

^ 126.0 

14.5 

27.0 

"31.7 

77.0 

23.0  r 

127.0 

15.1 

28.0 

26.3 

70.^ 

— 24.3/ 

128.0 

15.4 

20.0 

11.6 

70.0 

24.7' 

129.0 

15.4 

30.0 

14.0 

80.0 

' 25.0.'.  ' 

130.0 

15.4 

31.0 

. 17.6 

81.0 

25.8  - . 

131.0 

15.2 

32.0 

' 8,  a., 

82.0 

25.5 

•.  . / 132.0 

15.4 

33.0,/ 

2l;3 

'83.0 

25.7 

133.0 

15.1 

34.0 

42.4 

**  • 

— - 

84.0 

25.6 

134.0 

15.0 

35.0 

41.0 

85.0 

25.3 

135.0 

14.6 

36.0 

41.6 

86.0 

25.1 

136.0 

14.3 

37.0 

41.1 

87.0 

24.4 

137.0 

14.0 

38.0 

40.6 

88.0 

24.0 

138.0 

13.3 

30.0 

40.2 

80.0 

23.3 

139.0 

13.2 

40.0 

30.3 

00.0 

23.7 

140.0 

12.5 

41.0 

30.2 

01.0 

23.4 

161.0 

12.2 

42.0 

37.0 

92.0 

22.2 

142.0 

11.7 

43.0 

36.4 

93.0 

21.2 

143.0 

10.7 

44.0 

37.1 

04.0 

20.0 

144.0 

10.0 

45.0 

33.4 

95.0 

18.5 

145.0 

8.5 

46.0 

34.2 

06.0 

17.3 

146.0 

8.6 

47.0 

34.4 

07.0 

15.8 

147.0 

7.0 

48.0 

32.4 

98.0 

13.8 

148.0 

5.9 

40.0 

32.4 

90.0 

15.6 

149.0 

5.7 
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SELF -SCREENING  JAMMING  EXAMPLE 


APPENDIX  C 

SELF-SCREENING  JAMMING  EXAMPLE 
SELF -SCREENING  JAMMING  EXAMPLE 

In  this  example,  the  characteristics  of  the  radar,  threat  missile,  and  environment 
are  identical  to  the  scenario  in  appendix  B.  There  is  only  one  enemy  aircraft  that 
launches  the  threat  missile.  There  are  no  standoff  jamming  aircraft.  The  threat 
missile  has  a self-screening  noise  jammer  on  board  with  10  watts  of  power  over  a 
10  megahertz  bandwidth.  The  threat  missile  flies  the  same  flight  profile  as  in 
appendix  B. 

The  inputs  to  the  detection  model  are  presented  in  table  C-1.  Figure  C-1  pro- 
vides the  input  deck.  The  numerical  outputs  (table  C-2)  and  the  graphical  outputs 
(figure  C-2)  present  the  results  of  the  self-screening  jamming  example. 


16  2 15000  0 150  0 1 1000  35  25  1000  10  5 50  1 0.5  5 100  1 1 


1 1 1 -1.55605  0.54851  1 1 1 150  1 00000001  1.333 


SELF-SCREENING  JAMMER  CASE 


TABLE  C-1 

RADAR  DETECTION  MODEL  INPUTS 


IMPOT  DEFINITION  VALUE 

X(l)  aultipath  effect  Included  (l>yee;  0>no)  1 

X(2)  surface  clutter  Included  (i-yest  0«no)  1 

e 

X(3)  surface  clutter  coefficient  in  decibels.  -30 

A ratio  of  1 square  meter  of  clutter  to 
1,000  square  meters  illuminated  would  be 
entered  as  -30  (see  reference  1,  pages 
527-534) 


X(4)  weather  clutter  included  (l>yes;  0«no)  1 

X(5)  rainfall  rate  in  millimeters  per  hour  1 


X(6)  attenuation  included  (l>yes;  0»no)  1 

X(7)  and  X(8)  X(7)-L2^  and  X(8)<«L2  of  natural  logarithmic  -1.55605 

regression  L^L^+L-LnR,  where  L is  the 
attenuation  loss  in  decibels  and  R is  0.54851 

target  range  in  nautical  miles  (see  section 
on  attenuation 


X(9) 

X(10) 

X(ll) 

X(12) 

X(13) 

X(14) 


jamming  included  (l>yes;  0*no)  1 

type  of  jammer  (0>standoff  jammer;  1 

l>self-screening  jammer) 

location  of  jammming  (0-side-lobe;  1 

1-maln-lobe;  2-both  main-  and  side-lobes) 

range  from  radar  to  jammers  in  nautical  150 

sdles.  Note:  all  jammers  must  be  at 
the  same  range 


noise  jammer  bandwidth  in  megahertz.  Note:  10 

all  jammers  must  have  identical  bandwidths 

effective  radiated  power  of  the  jammer (s)  00000010 

in  watts.  Packed  word  "SSSSMMMM"  where 

"8SSS"  refers  to  side-lobe  jamming  and 

"MMMM"  refers  to  main-lobe  jamming.  Note: 

all  jamming  power  in  the  appropriate  lobe 

snist  be  totaled  and  entered  as  one  value. 

For  example,  three  100  watt  jammers  in  the 
side-lobe  and  two  100  watt  jammers  in  the  main- 
lobe  would  be  entered  as  03000200 


TABLE  G-1  (Cont’d) 


INPOT 

DEFINITION 

VALUE 

• 

X(15) 

coefficient  of  refractlvity  (1.333  for  the 

4/3  earth  approximation) 

1.333 

X(16) 

•ea  state  (Beaufort  scale) 

2 

X(17) 

tarset  altitude  at  the  beginning  of  the 
run  (in  feet) 

15,000 

» 

x(ie)  . 

target  altitude  at  the  end  of  the  run  (in 
feet) 

0 

X(19) 

target  range  at  the  beginning  of  the  run 
(in  nautical  miles) 

150 

X(20) 

target  range  at  the  end  of  the  run  (in 
nautical  miles) 

0 

§ 

X(21) 

average  target  cross  section  (in  square 
BMters) 

1 

X(22) 

Peak  radar  transmitted  power  (in  kilowatts) 

1000 

X(23) 

main-lobe  antenna  gain  (in  decibels) 

35 

s 

X(24) 

side- lobe  antenna  gain  (in  decibels) (e.g. , 
side-lobe  down  25  db  from  the  main-lobe 
gain  would  be  entered  as  25) 

25 

X(25) 

radar  frequency  (in  megahertz) 

1000 

X(26) 

receiver  noise  bandwidth  (in  megahertz) 

10 

X(27) 

receiver  noise  figure  (in  decibels) 

5 

X(28) 

integration  improvement  in  S/N  ratio  (in 
decibels) (e.g. , if  the  radar  integrates  100 
pulses,  perfect  integration  improvement  would 
be  entered  as  20) (see  reference  i.  pages 
35-40) 

50 

t 

X(29) 

radar  system  losses  (in  decibels) 

1 

X(30) 

radar  pulse  length  (in  microseconds) 

.5 

X(31) 

azimuth  beamwidth  in  degrees  measured  at 
the  3 db  down  level 

5 

1 

X(32) 

radar  antenna  height  (in  feet) 

100 

X(33) 

radar  polarization  (0-vertical;  1-horizontal) 

1 

X(34) 

output  in  graph  form  (1-yes;  0-no) 

1 

I 
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TABLE  G-2 

SELF -SCREENING  JAMMING  EXAMPLE 


0.  0 

61.0 

50.0 

32.1 

100.0 

14.3 

1.0 

71.1 

51.0 

31.4 

101.0 

15.6 

2.  0 

62.4 

52.0 

21.7 

102.0 

12.9 

0.0 

50.1 

53.0 

26.9 

103.0 

11.0 

4.0 

56.0 

54.0 

23.0 

104.0 

3.4 

5.0 

54.9 

55.0 

19.  3 

105.0 

”4.5 

r>.  0 

53.7 

56.  0 

23.0 

106.0 

”0.7 

7.0 

52.5 

57.0 

17.9 

107.0 

1.7 

Q.  0 

51.5 

50.0 

6.5 

108.0 

■ "12.2 

9.0 

9.3 

59.0 

16.6 

109.0 

"12.2 

1*0.0 

45.7 

60.0 

”9.6 

110.0 

”4.4 

11.0 

34.  7 

61.0 

”27.5 

111.0 

"15.3 

12.0 

30.0 

62.0 

”4.6 

112.0 

”7.1 

13.0 

14.4 

63.0 

”3.1 

113.0 

6.1 

14.0 

”1.9 

^ 

64.0 

"36.5  -f ; 

114;0” 

6.3 

15.0 

22.6  , 

"os.o 

,6.4: 

115.0  ' 

7.4 

16.0 

31.1  / / 

66\0 

11.0  

11 6.-0' 

9.0 

17.0 

“4.2  /, 

67.'0 

12.7 

117-.0 

15.1 

in.o 

13.1' 

68.0 

12.9 

110.0 

15.3 

19.0 

10.9 

69.0 

20.6 

119.0 

16.5 

27.0 

■’35.3 

.-7 

70.0 

-- 16.2  y 

120.0 

19.5 

21.0 

19.0 

71.0 

10.0 

121.0 

20.0 

22.0 

”13.0 

72.0 

24.6 

122.0 

20.5 

23.0 

20.2 

73.0 

24.5 

123.0 

22.2 

24.0 

"39.0 

/ 

74.0  V 

25.2 

124.0 

22.  0 

25.  0 

19.9 

75.0 

27.5 

125.0 

22.7 

26.0 

5.7 

76.0 

■ 24.7' 

126.0 

22.8 

27.0 

"37.7 

77.0 

■ 27.4  , 

127.0 

23.5 

23.0 

21.1 

78.0  - 

' . 27.9 

120.0 

23.9 

29.0 

6.  2 

79.  0 

20.4,  , '•  ' 

129.0 

23.9 

30.0 

3.0 

80.0  .. 

29.C  i ; 

130.0 

24.1 

31.0 

12.7 

01.0 

29.8 

131.0 

23.9 

32.0 

4.0,.  

82.0 

29.6 

132.0 

24.2 

33.0  / 

16.9  : 

83.0 

29.9 

133.0 

24.0 

34.0  ' 

30.3  ■ 

■ - 

84.0 

29.9 

134.0 

23.9 

35.0 

37.1 

85.0 

29.7 

135.0 

23.6 

36.0 

30.0 

86.0 

29.6 

136.0 

23.4 

37.0 

37.7 

87.0 

29.1 

137.0 

23.2 

30. 0 

37.5 

80.0 

20.7 

138.0 

22.5 

39.  0 

37.2 

89.0 

20.1 

139.0 

22.5 

40.0 

36.6 

90.0 

20.7 

140.0 

21.9 

41.0 

36.7 

91.0 

20.5 

141.0 

21.7 

42.0 

35.7 

92.0 

27.4 

142.0 

21.2 

43.0 

34.  3 

93.0 

26.5 

143.0 

20.3 

44.0 

35.2 

94.0 

25.5 

144.0 

19.7 

45.0 

31.0 

95.0 

24.0 

145.0 

10.2 

46.0 

32.0 

96.0 

22.9 

146.0 

18.4 

47.0 

33.  2 

97.0 

21.5 

147.0 

16.9 

40.0 

31.3 

90.0 

19.7 

140.0 

15.0 

49.0 

31.6 

99.0 

21.6 

149.0 

15.7 
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PI^OGRAM  OFTECT 


RAOAR  DETECTION  MODEL  INVOLVING  MULTIPATH  EFFECTS . JAMMING, SURF ACE 
CLUTTER,  AND  HEATHER  CLUTTER.  TH^  JAMMING  MAV  BE  FROM  A STAND-OFF 
JAMMER  OR  FROM  THE  TtRGET  I TSELF ( SELF-SCREENING  JAMMEPI.  THE 
TARGET  MAY  BE  INPUT  AT  ANY  ALT  IT  UOE/RANGE  AND  HAV*^  ANY  FLIGHT 
PROFILE  ON  A RADIAL  FROM  THE  DEFENDING  SHI^. 


RAOAR  DETECTION  MODEL  INPUT  VARIABLES 

*11)  MULTIPATH  effect  I NCLUOED ( 1 =YES , 0=NO) 

X«2I  SURFACE  CLUTTER  INCLUOEOJ 1=YES,  0=NO) 

X<3)  SURFACE  CLUTTEF  COEFFICIE NT (OB)  (EG.,  A RATIO  OF  1 SQUARE 
METE?  CLUTTER  TO  1000  SQUARE  METER  ILLUMINATED  MOULD  BE 
INPUT  AS  -30) 

*((,)  HEATHER  CLUTTER  INCLU0ED(  1=  YES,  0=N0) 

X(5)  RAINFALL  PATE (MM/HP) 

X(6)  ATTENUATION  INCLUDED ( 1= YE S,  0=NO) 

X(7)  LOSS(l)  (LOSS  = LOSS(l)  ♦ L0SS(2)  X LN  R 

SEE  RESEARCH  CONTRIBUTION  FOR  CALCULATION  OF 
LOSS(l)  AND  LOSS (2)  NATURAL  LOGARITHMIC 
X(8)  L0SS(2)  REGRESSION  COEFFICIENTS) 

X(9)  JAMMING  INCLUOEO(l=YES,  0=NO) 

X(10)  TYPE  OF  JAMMER(0»STANO-OFF  JAMMER,  laSELF-SCREENING  JAMMER) 
X(ll)  LOCATION  OF  JAMMING ( 0=SIOE-L03E , l=MAIN-LOBE,  2=BOTH) 

X(12)  RANGE  TO  JAMMER(NH) 

X(13)  NOISE  JAMMER  BA  NOW lOTH (MHZ) 

XdA)  EFFECTIVE  RADIATED  POMER  OF  JAMMER(H) 

(PACKED  MORD  SSSSMMMM  WHERE  SSSS  REFERS  TO  SIDE-LOBE 
JAMMING  AND  MMMM  REFERS  TO  MAIN-LOBE  JAMMING) 

X(15)  COEFFICIENY  OF  REFRACTION (1. 3333  FOR  4X3  EARTH  ASSUMPTION) 
X(16)  SEA  STATE (BEAUFORT  SCALE) 

X(17)  TARGET  ALTITUDE  AT  THE  BEGINNING  OF  THE  RUN(FT( 

X(18)  TARGET  ALTITUDE  AT  THE  END  OF  THE  RUN(FT) 

X(19)  TARGET  RANGE  AT  THE  BEGINNING  OF  THE  RUN(NM) 

X(20)  TARGET  RANGE  AT  THE  END  OF  THE  RUN(NM) 

X(21)  AVERAGE  TARGET  CRO SStfSECT ION (SQ . M.) 

X(22)  PEAK  TRANSMITTED  POHER(KH) 

X(23)  ANTENNA  GAIN,  MAIN-LOBE (OB) 

X(24)  ANTENNA  GAIN,  SIDE-LOBE (OB)  (EG.,  SIDE-LOBE  DOWN  2S0B  MOULD 
BE  INPUT  AS  25) 

X(25)  OPERATING  FREQUENC Y (MHZ) 

X(26)  RECEIVER  NOISE  BANOHIOTH ( MHZ ) 

X(27)  RECEIVER  NOISE  FIGURE(OB) 

X(28>  INTEGRATION  IMPROVEMENT  IN  S/N  RATI0(D3) 

X(29)  SYSTEM  LOSSES(OB) 

X(30)  PULSE  LENGTH(MICROSECONOS) 

X(31)  AZIMUTH  3EAMHI0TH(0CG> (AT  30B  DOWN  LEVEL) 

X(32)  ANTENNA  HEIGHT(FT) 

X(33)  POLAPIZATION(0=VERTICAL,  l=HORIZONTAL) 

X(34)  PLOTTER  OUTPUTdsYES,  0=NO) 
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C — 

COMMON  X (351  ,RNM  (350)  ,SGNSOa(350)  , TEMP  ( 10  ) . •JELT  AR  , 
l^lt*2.  ANGLEl.  ANGLE2tPhI,PHl0,PHH/,P.o.0,NlM0,RV.0.IST0P.PI 

C“  j 

C REAO  ALPHAMERIC  SUBJECT  LINF  FOP.  THE  OUTPUT  PLOT 

C- 


1 READ  10,  (TEHP(l),I=l,ia) 

10  FOPMAT  (IQAai 

IF  (EOF, 60)  100,20 
20  PRINT  1000, (TEMP(I) ,1=1,10 
1000  FORMAT  (1H1,10A5) 

C 

C READ  PARAMETERS  INTO  THE  PROGRAM 

C — 

CALL  nATA2(X,60, INO) 

IF  (INO.GT.O)  go  to  100 
PRINT  30,  ((I,X(I)),I=1,35) 

30  FORMAT  (5X,7 (I4*=*F10. A), I4*=*,F10. 7/5X,5 (IA*=*F10.4) ,I4*=*F10 , 
-2(l4*=*F10.4)/2(5X,d(I4*=*Fl0. 4) /) ,5X, 2 ( I4*=*Fio. 4) , I4*=*F15 . 12/) 
C 

C COMPUTE  SIGNAL-TO-NOISE  RATIO  AS  A FUNCTION  OF  RANGE 

C 

CALL  R A0RNG( ISTART) 

IF  (X(34)  .EO.O.)  GO  TC  1999 

C 

C DETERMINE  SCALING  FACTORS  FOR  THE  OUTPUT  PLOT 

C 

SNMAX  = SNHIN  = SGNSDB ( ISTOP) 

I9EGIN  = ISTART+l 

lENO  = ISTOP 

00  50  JJ=IBEGIN,IEND 

IF  (SGNSOa(JJ) .LT.SGNS0B( JJ'l) ) GO  TO  40 
IF  (SGNSOB(JJ).GT.  SNMAX)  SNMAX=SGNS09 (J J) 

GO  TO  50 

40  IF  (SGNSD9(JJ).GT.SGNS0B(JJ-1))  GO  TO  50 
IF  (SGNSOBC JJ) .LT.  SNMIN)  SNMIN=SGNS0B ( JJ) 

50  CONTINUE 

SNMIN  =((INTF(  SNMIN/10.))*10)-10. 

60  SGNSOBIISTOP  ♦!)=  SNMIN 

SNMAX  =((INTF(((  SNMAX-  SNMIN) /7 ) / 10 ) ) *10) ♦! 0. 

SGNSOBdSTOP  ♦■2)=  SNMAX 
RNMdSTOP  ♦DsO. 

IF  (ISTOP  .LE.50)  RNMdSTOP  +2)=5. 

IF  (ISTOP  .LE.50)  GO  TO  70 
IF  (ISTOP  .LE.lOO)  RNMdSTOP  ♦2)  = 10. 

IF  (ISTOP  .LE.lOO)  GO  TO  70 
IF  (ISTOP  .LE.150)  RNMdSTOP  ♦•2)  = 15. 

IF  (ISTOP  .LE.150)  GO  TO  70 
IF  (ISTOP  .LE.200)  RNMdSTOP  ♦2)=20. 

IF  (ISTOP  .LE.200)  GO  TO  70 
IF  (ISTOP  .LE.250)  RNMdSTOP  ♦2)=25. 

IF  (ISTOP  .LE.250)  GO  TO  70 
IF  (ISTOP  .LE.300)  RNMdSTOP  ♦2)=30. 

70  IPRINT  = IENO+2 

NUMBER* (ISTOP-ISTART)*! 
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DETECT  mns/iT 

INUM  » 1 

00  9999  JI  » ISTAPT.IPPINT 
9M>HIKUHl  » 9NM(JII 

SGNSOSdNUNI  3 SGNSOnt  Jll 
INUn  = INUH*! 

9999  CONTIHUE 

: PLOT  SIGN4L-T0-N0ISE  FATIO  AS  A FUNCTION  OF  RAN^E 

CALL  PL0TTEF.2  (RNN,  SGNSOB  , NUMBER  t 60 1 lAHRANGE  (N. 

ISHS/N  <0B),6tTEMP,S0,e) 

1999  CONTINUE 
PRINT  2000 

2000  FORMAT  < IIX, •PNG*, 7X ,*S/N* I 
00  SOOO  I s 1, NUMBER 

PRINT  3000,I,FNM(II ,SGNSOB(I> 

3000  FORMAT  ( I4,F 1 0 . 1 , FlO . l I 
9000  CONTINUE 
GO  TO  1 
100  END 
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SUBROUTINE  RAC<?Nr,(I3TfiRT> 


SUBROUTINE  TO  OETERHINF.  S IGNAL-T O-NOISE  RATIO  AS  A FUNCTION  OF  RANGE  | 

'|:i 

COMMON  X»35),  = Ni(350»«SGN'SDfl(3';0»tTEMP«10),OELTAR,  ' | 

IRltRZ, ANGLEl , ANGLE2 , PHI .PHIO.PHI V , RROt RHO.RV .D t ISTOP»PI  | 

DIMENSION  ALPHA1(2)  | 

OB»ZZ»=10*ALOG10I7ZI  | 

UNOBIZZZ)=10.**(ZZZ/1C) 


CONVERT  INPUT  VARIABLE  TO  UNITS  OF  MEASURE  NECESSARY  FOR  COMPUTATION 


XMPATM=X( 1» 

SURFCF=X(2) 

SIGO=UNOe(X<3ll 

HEATHR=X(4) 

RAIK=X(5) 

ATTNTN=Xt6) 

ALPHAKH  = X(7) 

ALPHAKZI  = X(8) 

XJAMsXI9> 

TYPJAM=Xfl<H 

STBJAM=Xtll» 

RJ=X(12»*1852.0035 

9J=X(13)*(10.**6» 

PJGJ=X(14» 

RE=20891199*XU5» 

NAVETM=10.*( (X(16»/10. » ••?) 

ALTBGN  = X«17» 

ALTEND=X(18» 

RNGBGN  = X(19»*6076. 11549 
RNGENO  = X«20»*6076. 11549 
SIGT=X(2U 
PT=XI22»*1000. 

GT=LN0B(X(23I ) 

GT1=UN0B(X(24)) 

10  GT1=GT/GTJ 
20  NAVELT=299.7925/X«25) 

9N=X(26»*(10.**6) 

FN=UN0B<X(27n 
XNElNsUNOSIX {28> ) 

SL=UNOB<X(29)) 

TAO=X<30)/tlO.**6» 

BZ=X(31»*. 017453293 
HlsX(32> 

RPOLAR=X<33» 

S1GJAM=0. 

OBOPsO. 

CLT=2.997925*«ia.**8» 

C=l. 38/110. ••23» 

TO- 290 

PI*3. 14159265 
AMBOA>HAVELT/. 3048006096 



C CALCULATE  PERMITTIVITY  AND  CONOUCTIVITY  CONSTANTS 


ooo  ooo  ooo  ooo 
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El  s 62.1521«583*2J.5€35195(»*AH0OA 

sigma  = 21.0811l081-<5  7.733»mi*AM30fl»*<<*7,8i338932*AM80A**2) 
FSIG»IPT*<GT**2t •XNEINI/C  n4*PI» ••3)*SLI 
FTAR*MAVELT»*2  . 


CALCULATE  THE  RAOAM  RECEIVER  NOISE 


RN0ISE=C*T0*FF«8N 
OENOM=RNOISE 
PkEALT  = ALT  END 

9ETA1  » IPt/EJ-fPNGOGN/JRE^ALTSGNt » 

0ETA2  = <pI/2)-(RNGEN0/(RE»ALT£N0) » 

GAMMA  = (ALT9&N-ALTEN0) /(BETA1-9ETA2» 

STEP  * (9ETA2-BETA1»/(X(19I-X  (20n 
I3TOP  » X(19» 

START  = (IALTEN3-Hl»/6076.1l549»»t.5 
IF  <X(20I  .GT.  STARTJ  ISTART  s X(20» 

IF  (X(20»  .LE.  START)  ISTART  = START 

DO  300  I = ISTART. I3TOP 

H2  = RE4PREALT 

SlOei  » (S1NF(BETA2) )*H2 

SIOEl  » SIOEl-(RE^Hl) 

SIOE2  * ICOSF(BETA2))*h2 

RT  = ((SIOE1**2I»(SIOE2**2))**Q.5 

RHORIZ  = 1.068*(SQRTF<X«15)l)*ISQRTF«Hl)*S(lRTF<H2-RE)l*b076.11549 
IF  CRT  .GT,  RHORIZ)  GO  TO  400 
R s RT»Q, 3048008095 
RR  : R/1852. 0035 

ALPHA  » ALPHAl(t)4’(ALPHAll2)»L0GF(RR)) 

ALPHA  3 ALPHA/ (RR«1852. 0035) 

ALPHA  * CALPHA/4.34)*(-1I 


CALCULATE  THE  SIGNAL  RECEIVED  FROM  THE  TARGET 


SIGTAR»FSIG*FTAR*SIGT/CR**4) 

IF  CATTNTN.EO.O)  go  to  30 
XNUM  3 SIGTAR*EXPFC2.*ALPHA*R) 
GO  TO  40 
30  XNUM3SIGTAR 


CALCULATE  THE  MULTIPATH  EFFECT 


40  IF  (KMPATH  .EO.  0.)  GO  TO  SO 

CALL  MPATH(RE.H1,H2,RT,AMB0A«E1,SI6HA,H,RP0LAR,0S0PI 


CALCULATE  THE  JAMMING  SIGNAL  STRENGTH 


50  IF  (XJAM.EQ.O)  GO  TO  80 

IF  CTTPJAM  ,LT,  .5)  GC  TO  60 

RJ  3 R 

60  IF  (STBJAH  .LT.  .5)  CO  TO  61 
IF  CSTBJAH  .LT.  1.5)  60  TC  52 
IF  (STBJAM  .LT.  2.5)  GO  TC  63 

61  CALL  SIDE  (P JG J, 6T1, HA VEL T ,BN,PI ,B J.R J. SIG JAM) 

e 
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GO  TO  70 

GZ  call  main  (PJGJtGT.HAVELT,QN,PI.BJ.RJ,SIGJAM) 

GO  TO  70 

61  CALL  BOTH  (P JGJ, GT , GTl « HA VELT , BN ,PI , BJ.RJ. 3IGJAHt 

70  IF  (ATTNTN  .EO,  01  GO  TO  71 
SI6JAM  = SIGJAM*EXPF(ALPHA*RJ» 

71  OENOM  = OENOMfSIGJAH 

C 

C CALCULATE  THE  SURFACE  CLUTTER  EFFECT 

C 

80  IF  (SURFCE.EQ. 0)  GO  TO  120 

XRNG=5639.039«*(SORTF(X«16m*(SQRTF(Hl>) 

IF  (RT.LE.XRNGI  GO  TO  90 
FACT2=0. 

GO  TO  110 

90  FACT2=IPT*<GT**2»*(WAVELT)*8Z*CLT*TAU*SIG0»/ 

1 <2. *(<«». •PT»**3»  I 
IF  (ATTNTN. EQ. 01  GO  TO  100 
FACT2  = (FACT2*EXPF(2.*ALPHA*9H /4R**3) 

GO  TO  110 

100  FACT2=FACT2/ (P**3) 

110  0EN0H=DEN0M»FACT2 

c 

C CALCULATE  THE  HEATHER  CLUTTER  EFFECT 

C 

120  IF  (HEATHR.EQ.OI  GO  TO  1L5 

FACT3=(186.*PT*GT*CLT*TAU*<PI • (RAIN**1. 6H / 

1(128.*  (MAI/ETH**2»*(10.  **181  » 

IF  (ATTNTN. EQ.OJ  GO  TO  130 
FACT3=FACT3*EXPF(2.*ALPHA*R)/(R**2» 

GO  TO  l<ia 

130  FACT3  = FACT3/(R**2» 

140  OENOH  = 0EN0M>FACT3 

149  SGNS  = XNUH/OFNOH 
SGNS-OB(SGNSI 
OBDP  = 0B0P*2. 

190  SGNSDa(II=3GNS*0B0P 
RNHdl  = P/1852.  0035 
GO  TO  200 

9000  PRINT  9010«RNM(I» 

9010  FORMAT  (5X,*TARGET  RAKGE=* ,E17.9» 

9011  PRINT  9020.SGNSOB(II 

9020  FORMAT  (5X,*TOTAL  S/N  RAT I0=*.E1 7. 91 

9021  IF  (XMPATH  .EO.  0.)  GO  TO  9061 

9022  PRINT  9030, SGNS, 090P 

9030  FORMAT  (5X,*S/N  PATIO  H/0  MULTIPATH**, El7, 9, 5X, *MULTIPATM  EFFECT** 
1,E17.9) 

9031  PRINT  9040,RT,R1,R2,OELTAP,ANGLE1,ANGLE2 

9040  FORMAT  (5X , *RT=* ,F10 . 2, 5X , *R1=*, FIO. 2, 5X , *R2»*,F10 . 2,5X , *OELTAR**, 
1F10.2,5X,*ANGLE1=*,E17.9,5X,*ANGLE2**,E17.9» 

9042  PRINT  9090,RRO,RHO,RV,O 

9090  FORMAT  (9X,*MULTIPATH  C0EFF=*,E17. 9,9X,*R0U6H  SEA  COEFF**, E17. 9/, 9 
1X,*SM00TH  SEA  C0EFF=*,E17,9,5X,*0IVERGENCE**,E17.9» 

9092  PRINT  9060, PHI, PHI0,PHIV 

9060  FORMAT  (5X,*T0TAL  PHASE  OIFFs*,E 17.9, 9X, •RANGE  PHASE  0IFF**,E17.9, 
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15it,*REFUECTE0  PMAGE  01  FF»  • ,E  17. 

GO  TO  9063 

9061  PRINT  9062 

9062  FORMAT  <5X,*MULT IPATH  EFFECT  NOT  CONSIOEKED’A 

9063  IF  (XJAM  .EO.  0.1  GO  TO  9092 
906U  IF  <ATTNTN  .EO.  1.)  GC  TO  906? 

9065  PRINT  9080,SIGJAM 

9060  FORMAT  (5V.«UNATTENUATE0  JAMMING  SIGNAL**.E17.9I 

9081  GO  TO  9101 

9082  PRINT  9090.SIGJAM 

9090  FORMAT  (5 X, • AT TENUA TED  JAMMING  SIGNAt»*.E17. 9» 

9091  GO  TO  9101 

9092  PRINT  9100 

9100  FORMAT  (5X,*NC  JAMMING  SIGNAL*! 

9101  IF  JSURFCE  .EO.  0.)  GO  TO  9122 

9102  IF  (ATTNTN  .EO.  1.1  GC  TO  9112 

9103  PRINT  9110, FACT? 

9110  FORMAT  «5X,»U!. ATTENUATED  SURFACE  CLUTTER»*,E17. 91 

9111  GO  TO  9131 

9112  PRINT  9120, FACT? 

9120  FORMAT  <5X ,* ATTENUATED  SURFACE  CLUTTERa*,E17.9! 

9121  GO  TO  9131 

9122  PRINT  9130 

9130  FORMAT  J5X, •SURFACE  CLUTTER  NOT  CONSIOEPEO*! 

9131  IF  (MEATHR  .EO.  0.1  GC  TO  9152 

9132  IF  UTTNTN  .EO.  1.)  GO  TO  9142 

9133  PRINT  9140, facts 

9140  FORMAT  (5X,*UNATTENUATEO  WEATHER  CLUTTER**, E17. 9// I 

9141  GO  TO  200 

9142  print  9150, facts 

9150  FORMAT  (5X, •ATTENUATED  WEATHER  CLUTTER, •,E17 . 9//I 

9151  GO  TO  200 

9152  PRINT  9160 

9160  FORMAT  (5X, •HEATHER  CLUTTER  NOT  CONSIOEREO^F/I 
200  9ETA2  * BETA2-STEP 

PREALT  * PREALT- (GAHMA^STEPI 
OENOH  s RNOISE 
300  CONTINUE 
GO  TO  600 
400  PRINT  500 

500  FORMAT  (5X, •TARGET  9EVONO  RADAR  HORIZON*! 

600  RETURN 
ENO 
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SUBPOUTINF  MPtTM(Re»Hl,M2,RT,AMB04fEltSIGMft,H,RPOLARt090PI 


SUBROUTINE  TO  DETERMINE  THE  MULTIPATH  EEFEDT 


COMMON  X(35) ,RNM (350 ) . SGNSOB (350  I . TEMP ( 10) . OELTAR, 
lRl.k2, ANGLEl,AN6LE2,PHI«PHl0,PHIV,RRD,RH0,Rtf .O.ISTOPtPI 
COMPLEX  CFNSQO.CFN.E.CFFCNT.BE 
HT  = H2-RE 

FNUM=( (RE*Hl)**2l+( (RE+HT )**2)“(RT**2) 

FDEN0M=2* (RE*H1) •(RE»HT) 

ANS=FNUM/FOENCM 
THETA=AC0SF( ANS) 

XLEFT=0. 

RIGHT=THET4 

THETA1=THET4/2.  ^ 


ITERATIVE  PROCESS  TO  CALCULATE  ANGLE  OF  INCIOANCE/REFECTION 


DO  30  I=l»31  □ 

THETA2=THETA-THETA1 

CALL  DIST(Hl,THETAl,Rf ,R1) 

CALL  0IST(HT,THETA2,RE,R2) 

RSULT1=( (H1**2)^(2.*RE*H1)-(R1**2) )/(2.*RE*Rl) 

PSI1=ASINF(RSULT1)  $ ANGLEl=PSI 1* ( 360./ (2. •PI I ) 

RSULT2= ( (HT**2)* (2.*PE*HT)-(R2**2) )/(2.*RF*R2l 
PSI2=ASINF(RSULT2)  $ ANGLE2=PSI2* ( 360. / (2. •»! I ) 

IF  (PSI1.GT.PSI2)  GO  TO  10 
RIGHT=THE T41 
GO  TO  20 
10  XLEFT=THETA1 

20  THETA1=XLEFT*( (RIGHT-XLEFT)/2.I 
30  CONTINUE 


CALCULATE  PATH  LENGTH  DIFFERENCE 


DFLTAR=(R1»R2(-RT 


CALCULATE  PHASE  SHIFT  DUE  TO  PATH  LENGTH  DIFFERENCE 


PHI0*(DELTAR*2.*PI)/AMB0A 


CALCULATE  SIGNAL  STRENGTH  REDUCTION  AND  PHASE  SHIFT  CAUSED  BY 
REFLECTION 


G=((PI*H*SINF(PSIl)/AMaOA)**2) 

RH0=EXPF(-8*GI 

CFNSaO>(El*( 1. .0. I >- ( ( 0., 1. (•IS. 3*SIGHA»AM90A) 
CFNzCSQRKCFNSQO) 

E>SINF(PS11I«(1. .0.) 

IF  (RPOLAP.EO.OI  GO  TO  50 
CFFCNTx(E-CFN) /(E»CFN) 

GO  TO  50 


40  aEx(1..0.)/CFN 

CFFCNT*  ( E-8E  I / (E  ♦•)?) 
50  RV*CABS(CFFCNT) 
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PMIV«CANG(CFPCNT  > 

c - 

c calculate  total  phase  shift 

c 

PHI=PHI0-PH1 V 

/<PF*|P1*F  2»*F I 
0=1./SQRTF(1.*H» 

RR0=RH0*RV*0 



c calculate  hultioatm  effect 

c 

DELTAFs1.»(RPD**2»»<2.*RPC*':OSF(PHIM 

OBOP*10.*ALOG10<nELTAF» 

RETURN 

ENO 
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SUa^OUTIHE  OISTiXH.TM.RE, 


SUflROUTINF  TO  CALCULATE  DISTANCE  RAOA?  TO  POINT  OF  REFLECTION 

ANO  FRON  POINT  OF  REFLECTION  TO  TARGET 

IF  ITH.GT..0C7»  GO  TO  10 
0=.5*CTH**2» 

GO  TO  20 
10  0=1-C0SF(TH» 

20  T*RE*(RE+XHI 
T=T*2.*0 
T=T*IXH**?> 

XR»SQRTF<T» 

RETURN 

ENO 
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SU5R00riNE  HAIM  (0 JG J,GT« WAVCLT . AN,<>I ,0 J,R J,SIG JAHI 


CALCULATE  HAIL-L03E  JAHHIKG  SIGNAL 


IPJGJ  » PJCJ 

NPJGJ  > IPJGJ/10000 

HPJGJ  « IPJGJ-(10000«NPJGJI 

SIGJAM  = (HPJGJ*GT*(HAVELT««2)»BNI/(((A.*PI)»*2I«BJ»(RJ»«2I) 

RETURN 

END 


APPENDIX  E 


REFLECTED  RAY  PATH  LENGTH  CALCULATIONS 


tO 


In  the  section  on  multipath  effects,  it  was  claimed  that  the  technique  used  to 
determine  the  length  of  the  reflected  ray  path  calculated  the  point  of  incidence/re- 
flection to  within  one  inch  for  targets  with  a range  up  to  300  nautical  miles.  The 
portion  of  the  radar  detection  model  that  calculates  the  reflected  ray  path  length  was 
converted  to  the  APL  computer  language.  This  version  is  presented  in  aimex  E-1. 
Additional  programming  was  inserted  to  determine  the  error  in  the  point  of  reflection. 
Tables  E-1  through  £-4  present  the  results  of  various  target  altitudes,  ranges,  and 
flight  profiles.  The  error  between  the  calculated  point  of  reflection  and  the  actual 
point  of  reflection  was  determined  to  be  less  than  one  inch. 


TABLE  E-1 

REFLECTED  RAY  PATH  LENGTH  CALCULATIONS 


TARGRT  ALTITUDE 

AT 

THE  BEGIRRIRG  GPP  THE 

RUR  *300  PEET 

TARGET 

ALTITUDE 

AT 

THE  ERD  GPP 

THE  RUR 

• 300  PEET 

• TARGET  RARGE  AT 

THE 

' BEGIRRIRG 

OP  THE  RUR 

m 25  RAUTICAL 

MILES 

TARGET 

RARGE  AT 

THE 

' ERD  OP  THE  RUR 

> 0 RAUTICAL 

MILES 

RADAR  ARTERRA  HEIGHT 

• 60  PEET 

TARGET 

R\ 

J72 

DELTA 

ARGLE  DISTARCE 

RARGE 

RARGE 

DIPPERERCE 

ERROR 

(RM) 

{.FEET) 

iPSET) 

iPEET) 

{DEG) 

(IRCH) 

24.9 

41136.3 

110766.0 

8.5ff"2 

2.2026824807~12 

0._70 

23.9 

38486.8 

107339.3 

1.7S"2 

7.3896444517“13 

0._70 

22.9 

35920.2 

103829.9 

5.7ff"2 

1.3642420527"12 

0._70 

21.9 

33441.8 

100232.3 

1.3£"1 

1.9610979507“12 

0._70 

20.9 

31056.1 

96542.0 

1.6ff“l 

1.6484591467“12 

0._70 

19.9 

28766.8 

92755.2 

1.4£"1 

6.6222582977“12 

0._70 

18.9 

26576.8 

88669.1 

5.8«"2 

8.5265128297"14 

0._70 

17.9 

24487.5 

84882.3 

9.7F"2 

2.5579538487“12 

0._70 

16.9 

22499.4 

60794.5 

1.7ff“l 

1.7621459847"12 

0.„70 

15.9 

20611.3 

76606.5 

1.9fr"l 

1.7053025657"13 

0._70 

14.9 

18821.2 

72320.7 

3.6/f“l 

1.8189894037"12 

0._70 

13.9 

17125.5 

67940.2 

2.0E~1 

3.2969182937~12 

0._70 

13.0 

15520.0 

63469.7 

2.0E~1 

.2.6147972677"12 

0._70 

12.0 

13999.3 

58914.4 

2.SE~1 

3.0695446187"12 

0._70 

11.0 

12557.6 

54260.1 

3.4J?~1 

7.8443918027“12 

0._70 

10.0 

11188.5 

49573.5 

6.9P~1 

4.7748471847"12 

0._70 

9.0 

9885.3 

44800.7 

7.5S"1 

6.8212102637“12 

0._70 

8.0 

8641.1 

39968.8 

6.5R“1 

3.9790393207"12 

0._70 

7.0 

7449.1 

35085.0 

8.5«"l 

3.1832314567~12 

0._70 

6.0 

6302.3 

30155.8 

8.8ff"l 

9.7770680447“12 

0._70 

S.O 

5194.1 

25188.2 

1.0J?0 

2.2509993867“11 

0._70 

4.0 

4117.8 

20189.5 

1.7ff0 

6.59‘38365877~12 

0._70 

3.0 

3066.6 

15165.1 

1.8?0 

3.1377567217~11 

0._70 

2.0 

2034.6 

10123.2 

3.370 

1.7280399337“ll 

0._70 

1.0 

1015.5 

5071.4 

6.070 

2.0918378147~11 

2.37" 
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1 
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TABLE  E-2  j 

REFLECTED  RAY  PATH  LENGTH  CALCULATIONS  I 


TARGST 

ALTITUDE 

AT  TBE^BEGIRRIRG  OF  THE 

RUR  • 10,000  FEET 

TARGET 

ALTITUDE 

AT  THE  ERD 

OF 

THE  RUR 

■ 10,000  FEET 

• TARGET 

RARGE  AT 

THE  BEGINRIRG 

OF  THE  RUR  » 100  RAUTICAL  MILES 

TARGET  RABGE  AT  THE  END  OF 
RADAR  AETERRA  HEIGHT 

THE  RUR 

■ 75  RAUTICAL  MILES 

m 60  FEET 

TARGET 

RARGE 

R1 

R2 

DELTA 

RARGE 

ARGLE 

DIFFERERCE 

DISTARCE 

ERROR 

(RM)  i 

FEET) 

(FEET) 

(FEET) 

(DEG) 

(IRCR) 

99.9 

9720.5 

597852.6 

7.0ff"l 

8.6401996667*11. 

0._70 

98.9 

9333.7 

592165.6 

7.3ff”l 

5.7980287237“12 

0._70 

97.9 

8970.3 

586455.2 

7.6P”1 

6.7188921097"ll 

0._70 

96.9 

8628.3 

580723.4 

8.1S“l 

9.4473762147*11 

0._70 

9S.9 

8306.1 

574971.7 

8.2F"l 

4.7293724497*11 

0._70 

94.9 

8002.1 

569202.0 

8.5ff"l 

8.98126018.7*12 

0._70 

99.9 

7714.8 

563415.5 

8.8»“l 

1.0905978347*9 

0._70 

92.9 

7443.1 

557613.5 

9.3ff"l 

4.3883119367*11 

0._70 

91.9 

7185.6 

551797.3 

9.8ff“l 

1.5347723097*10 

0._70 

90.9 

6941.4 

545967.8 

l.OFO 

1.5143086787*10 

0._70 

69.9 

6709.4 

540126.0 

l.OFO 

1.2050804797*10 

0._70 

86.9 

6488.6 

534272.9 

1.070  . 

1.6234480427*10 

0._70 

87.9 

6278.9 

528409.2 

1.170 

1.5575096767*10 

0._70 

86.9 

6078.7 

522535.7 

1.170 

8.0490281107*11 

0._70 

85.9 

5887.8- 

516653.0 

1.170 

6.8212102637*12 

0._70 

84.9 

5705.4 

510761.8 

1.270 

1.7826096157*10 

0._70 

89.9 

5531.0 

504862.6 

1.270 

1.9144863477*10 

0..70 

83.0 

5364.1 

498955.9 

1.370 

9.0767571237*10 

0._70 

82.0 

5204.2 

493042.3 

1.370 

1.2187229007*10 

0._70 

81.0 

5050.8 

487122.1 

1.370 

7.9580786407*11 

0._70 

80.0 

4903.6 

481195.8 

1.470 

9.9589669847*11 

0._70 

79.0 

4762.2 

475263.7 

1.470 

1.1914380597*10 

0._70 

78.0 

4626.2 

469326.3 

1.570 

2.1850610217*10 

0._70 

77.0 

4495.3 

463383.8 

1.570 

3.1673152987  10 

0._70 

76.0 

4369.2 

457436. 5 

1.670 

1.1777956387  10 

0..70 
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TABLE  E-3 

REFLECTED  RAY  PATH  LENGTH  CALCULATIONS 


TARGET  ALTITUDE  AT  THE  BEGIHHIHG  OF  THE  HUH  = 
TARGET  ALTITUDE  AT  THE  EHD  OF  THE  RUN  s 
TARGET  RANGE  AT  THE  BEGINNING  OF  THE  RUN  = 
TARGET  RANGE  AT  THE  END  OF  THE  RUN  * 
RADAR  ANTENNA  HEIGHT  « 


75,000  FEET 
75,000  FEET 

100  NAUTICAL  MILES 
75  NAUTICAL  MILES 
60  FEET 


TARGET 

R1 

R2 

DELTA 

ANGLE 

DISTANCE 

RANGE 

RANGE 

DIFFERENCE 

ERROR 

(.NM) 

(FEET) 

(FEET) 

(FEET) 

(DEG) 

(INCH) 

100.6 

535.7 

610871.7 

1.3?! 

1.121043169ff“8 

4.0i?“6 

99.6 

529.4 

604856.8 

1. 3ffl 

3.754394128/?“9 

1.3ff“6 

98.6 

523.2 

598842.7 

1.3jri 

1.424814399ff~8 

5.2£:"6 

97.6 

517.0 

592829.6 

1.3^1 

2.213710104i?“9 

8.3ff"7 

96.6 

510.8 

586817.4 

1.4ffl 

4. 549292498ff~9 

1.5ff"6 

95.6 

504.6 

580806. 2 

1.4^1 

7.870767,a49ff“9 

2.8F"6 

94.6 

498.5 

574796.0 

1.4i71 

3.170498530?~9 

1.1£"6 

93.6 

492.4 

568786.9 

1.4£:i 

3.153218131P"8 

l.lff*5 

92.6 

486.4 

562778.8 

1.4ffl 

1.006083039ff”8 

3. 5r"6 

91.7 

480.3 

556771.9 

1.4ffl 

1.256194082F"8 

4. Sff"6 

90.7 

474.3 

550766.2 

1.5ffl  . 

9.014911483?"9 

Z.2E~6 

89.7 

468.4 

544761.6 

1.5F1 

3.218519850ff~8 

1.1E~5 

88.7 

462.4 

538758.3 

1.5F1 

3.432433004ff“8 

1.2E~S 

87.7 

456.5 

532756.3 

i.5;ri 

2. 300293999?:"8 

9.2E~6 

86.7 

450.7 

526755.6 

1.5?1 

1.511580194K~8 

5.4ff~6 

85.7 

444.8 

520756.3 

1.6ffl 

7.777998689ff“9 

2.SE~6 

84.7 

439.0 

514758.4 

1.6i?l 

3.105014911ff“8 

1.1E~5 

83.8 

433.2 

508762.0 

1.6ffl 

7. 305061444P“9 

2,6E~6 

82.8 

427.4 

502767.2 

1.6P1 

3.  342938725ir“8 

I.IE'S 

81.8 

421.7 

496773.9 

1.7i71 

2.122033038r"8 

7,6E~6 

80.8 

416.0 

490782.3 

1.7F1 

1.455555320.r“8 

5.3E~6 

79,8 

410.3 

484792.4 

1.7E1 

2.656088327i:~8 

9.5E~6 

78.8 

404.6 

478804.4 

1.7ffl 

1.958323991F"8 

7.0E~6 

77.8 

399.0 

472818.1 

1.8ffl 

7.490598363P"9 

2 . 7,'7~6 

76.8 

393.4 

466833.8 

1.8^1 

4.152389010ff“8 

1.4r"5 

E-4 


6 


f I 

I 

» f 


TABLE  E-4 

REFLECTED  RAY  PATH  LENGTH  CALCULATIONS 


TARGET  ALTITUDE  AT  THE  BEGIHHIHG  OF  THE  RUH  = 75,000  FEET 

TARGET  ALTITUDE  AT  THE'  END  OF  THE  RUH  = 75,000  FEET 

TARGET  RANGE  AT  THE  BEGIHHIHG  OF  THE  RUH  ~ 25  NAUTICAL  MILES 

TARGET  RANGE  AT  THE  END  OF  THE  HUH  = 1 NAUTICAL  MILE 

RADAR  ANTENNA  HEIGHT  * 60  FEET 


I. 


ARGET 

J?1 

R7 

DELTA 

ANGLE 

DISTANCE 

ANGE 

RANGE 

DIFFERENCE 

ERROR 

(NM) 

(FEET) 

(FEET) 

(FEET) 

(DEC) 

(INCH) 

27.8 

136.0 

169116.7 

5.3E1 

1.356238499S"8 

4. 8E~8 

26.9 

131.6 

163704.6 

5.4?1 

3.508466761ff"8 

1.2J?“5 

26.0 

127.2 

158339.6 

5.6£ri 

5.4104020Elff~8 

1.9F“5 

25.1 

122.9 

153026.8 

5.8J71 

5.935726221ff~8 

2.ir“5 

24.3 

118.6 

147771.8 

6.o;ri 

5.321635399P~8 

1.9E~S 

23.4 

114.4 

142581.0 

6.2?1 

7.172639016ff“8 

2.5E"5 

22.6 

110.3 

137461.4 

6.5£1 

1. 308508217ff~7 

4.7ff"5 

21.8 

106.2 

132421.7 

6.7F1 

2.437445800B- 8 

8.7F"6 

20.9 

102.2 

127471.1 

7.0E1 

1.16779119'. "7 

4.2P"5 

20.1 

98.3 

122620.5 

7.3E1 

1.156149664S“7 

4.1ff"5 

19.4 

94.4 

117882.2 

7.6E1 

1.378066372i?“8 

4.9ff"6 

18.6 

90.7 

113270.2 

7.9E1 

7.43602868ir~8 

2*.6ff"5 

17.9 

87.1 

108800.7 

8.2E1 

• 8.221832104K“9 

2.9£:"6 

17.1 

83.6 

104491.9 

8.SE1 

6.303889676F"8 

2.2ff"5 

16.5 

80.3 

100364.6 

8.9E1 

2.908927854P"8 

1.0ff"5 

15.8 

77.2 

96442.1 

9.3«:i 

2.823071554ff"9 

1.0ff"6 

15.2 

74.2 

92750.4 

9.6E1 

1.290463842E~7 

4.6K"5 

14.6 

71.4 

89318.1 

1.0E2 

8.802453521ff"8 

3.1ff"5 

14.1 

68.9 

86176.1 

1.0E2 

1.05268S547;?"7 

3.7P”5 

13.7 

66.7 

83357.1 

1.0E2 

3.492459654P”8 

1.2ff"5 

13.3 

64.7 

80895.3 

1.1E2 

8.806819^)96^“8 

3.1F"5 

12.9 

63.0 

78823.7 

1.1E2 

8. 358620107ff“8 

3.0ff“5 

12.6 

61.7 

77174.0 

1.1E2 

6.443588063P"8 

2.3ff"5 

12.4 

60.7 

75974.0 

1.1E2 

1.725857146ff~8 

6.2ff"6 

12.3 

60.1 

75244.8 

1.1E2 

1.030275598ff”8 

3.6ff"6 

\ tC 


ANNEX  E-1 

REFLECTED  RAY  PATH  LENGTH  PROGRAM 


f 

tjtnon 

Cl] 

* TARGET  R\  E7 

DELTA 

AWGLE 

C2] 

*RARGE 

RAEGE 

DIFFEREECE 

[a] 

• (Wf)  (.FEET)  (FEET) 

(FEET) 

(DEG) 

Co] 

PI«-3.101S926S 

Cs] 

M«-20891199>i210 

C6J 

ALTBGE*Xi^ 

C7] 

ALTERD*Xn 

C8] 

1 8 X 6 0 7 6 . 1 1 S 0 9 

C9]. 

P 1 9 « 6 0 7 6 . 1 1 5 0 9 

CIO] 

ffl»231 

Cll] 

FREALT*AITBGR 

C12) 

BETA1*(PI*2)-(RRGBCR*(ER*ALTBGR)) 

Cl3] 

BETA 2*(PI*2)- (RRGEND*( ER*ALTEHD ) ) 

CIO] 

GAM14A*(AL TBGR -AL TERO ) *(BBTA1-BETA2) 

CIS] 

STEP*( BETA  2 -BETA 1)*(/18-219) 

CIS] 

AC-.H2*ERtPREALT 

C17] 

SIDEi*( lo(nETA 1))»R2 

CIS] 

SIDE1*SIDE1 -(ER*ni) 

C19] 

SIDE2*(20(BETAi))»H2 

C20] 

RT*((SIDEl*2)*(SIDE2*2))*0.S 

C21] 

RT-R2-ER 

C22] 

FRUH*((KR*R\)»2)*((ER*HT)*2)-(RT*2) 

[23] 

FDER0U*2 « ( ER*H1 ) * ( ER^RT) 

C20] 

AKS*PRUH*FDEROH 

C2S] 

TRETA*~2oARS 

C26] 

iEFT*0 

C27] 

RIGRT*TRE?A 

[28] 

VHSTAl*TasrA*2 

C29] 

I*i 

C30] 

BB : TRETA 2*?RETA - TRET A 1 

C31] 

DISTAHCE  TRETAl 

[32] 

R2*RT  DISTAHCE  TRETA2 

C33] 

RSUl.T\*(Rl»2)*(2»En»Rl)-(Ri*2) 

[30] 

RSVLT\*-RSVLT\  *(2i>ER*R\ ) 

[35] 

PSIi*~loRSULTl 

C36] 

RSUDT2*(HT*2)*(2»ER»RT)-(R2i>2) 

[37] 

RSULT2*-RSULT2  » ( 2 nr.R^R2  ) 

C38] 

PSI2*-~\oRSVLT2 

C39] 

*BCnxPSri>PSI2 

[00] 

RIGHT*TRETAl 

[OJ] 

*BD 

[02] 

BCi LEFT*TRETA\ 

[03] 

BDtTRETAl*LEFT*( (RIGRT-LEFT)*2 ) 

• 

[00] 

[OS] 

-*8B«i/s32 

[OS] 

DELTAR‘^(Ri.*R2)-RT 

[07] 

>( 1 « ( 36 0 ♦ ( 2 xPI ) ) 

[OS] 

AHGLE2*PSI2»(i(>0*(2»PI)) 

[09] 

DIFF*\ (PSI1-PSI2) 

[SO] 

*ZA»\PSI1>PSI2 

CSX] 

HEVPSZ\*PSI\*(DIFFt2) 

0Z8TAWCE 

mttojf 

uncH)* 


E-6 


*ZB 

ZAtllSUPSn*PSIl-(DIFF*2) 

ZB  tHISS*  I ( ( ( 20PSJ1 ) «;?!  » 1 2 ) - ( ( lOBSVPSiX  ) «iri « 1 2 ) ) 
Xl4-(/r7*6076.11S49)ff6  3 2 0 
X2«-i71ffl3  7 2 0 
X3-^ff2013  7 2 0 
K>i*I)BLTAB§\X  12  3 
Ki^DIFFtXZ  1 10  3 
Kt*MISSI>X2  2 2 3 
(X1).(X3).(X3).(X4).(XS).(X6) 

' PBBALT*PBBALT*{  GAMMAytSTEP) 

BETAX^BETAX*STEP 

EIID*{HX*Elf)*RT-H2 

•*0«iXX2)<0 

•*^0»\BKTAX*BSTA2 

*AC 

DSLTAH*iRX*R2)-RT 


0 R*FDR  DISTARCE  FDTR 

Cl]  ■*CXKtXDrX>0.007 

C2D  FDD-0. FDTR *2) 

1 33  -LB 

t4]  GRtFDD-X-{2oFDTH) 
t5]  LRiFDT-RR*{ER*FDn) 

[6]  FDT-FDT»2*FDD 

[7]  FDT*FDT*IFDR*2) 

Ca]  X*FOr*0,5 

7 


